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^ I In this paper, we propose a transmit beamspace energy focusing technique for multiple-input 

multiple-output (MIMO) radar with application to direction finding for multiple targets. The general 
angular directions of the targets are assumed to be located within a certain spatial sector. We focus 
the energy of multiple (two or more) transmitted orthogonal waveforms within that spatial sector using 



(N 



f/3 ' transmit beamformers which are designed to improve the signal-to-noise ratio (SNR) gain at each receive 

antenna. The subspace decomposition-based techniques such as MUSIC can then be used for direction 
finding for multiple targets. Moreover, the transmit beamformers can be designed so that matched- 
vQ ' filtering the received data to the waveforms yields multiple (two or more) data sets with rotational 

^^ \ two data sets or parallel factor analysis (PARAFAC) for more than two data sets. Unlike previously 



o 
o 



invariance property that allows applying search-free direction finding techniques such as ESPRIT for 



reported MIMO radar ESPRIT/PARAFAC-based direction finding techniques, our method achieves the 
rotational invariance property in a different manner combined also with the transmit energy focusing. 
As a result, it achieves better estimation performance at lower computational cost. Particularly, the 
proposed technique leads to lower Cramer-Rao bound than the existing techniques due to the transmit 



X 

^ . energy focusing capability. Simulation results also show the superiority of the proposed technique over 

the existing techniques. 
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I. Introduction 

The development of multiple-input multiple-output (MIMO) radar is recently the focus of 
intensive research [[I]|-|l3l. A MIMO radar is generally defined as a radar system with multiple 
transmit linearly independent waveforms and it enables joint processing of data received by 
multiple receive antennas. MIMO radar can be either equipped with widely separated antennas 
[l2l or colocated antennas Q. 

Estimating direction-of-arrivals (DOAs) of multiple targets from measurements corrupted by 
noise at the receiving array of antennas is one of the most important radar applications frequently 
encountered in practice. Many DOA estimation methods have been developed for traditional 
single-input multiple-output (SIMO) radar Hl- [[T2ll . Among these methods the estimation of 
signal parameters via rotational invariance techniques (ESPRIT) and multiple signal classification 
(MUSIC) are the most popular due to their simplicity and high-resolution capabilities O, [HI, 
[fT2l . Moreover, ESPRIT is a special and computationally efficient case of a more general 
decomposition technique of high-dimensional (higher than 2) data arrays known as parallel 
factor analysis (PARAFAC) ifBll. lfT4l. 

More recently, some algorithms have been developed for DOA estimation of multiple targets in 
the context of MIMO radar systems equipped with M colocated transmit antennas and N receive 
antennas [[T5l - [fT9l . The algorithms proposed in ifTSl and lfT6l require an exhaustive search over 
the unknown parameters and, therefore, mandate prohibitive computational cost if the search 
is performed over a fine grid. On the other hand, the search-free ESPRIT-based algorithms 
of [fTTl and [[TSl as well as PARAFAC-based algorithm of [19J utilize the rotational invariance 
property of the so-called extended virtual array to estimate the DOAs at a moderate computational 
cost. It is worth noting that in the case of MIMO radar, the advantages of the aforementioned 
DOA estimation methods over similar MUSIC- and ESPRIT/PARAFAC-based DOA estimation 
methods for SIMO radar appear due to the fact that the extended virtual array of MN virtual 
antennas can be obtained in the MIMO radar case by matched- filtering the data received by 
the A^-antenna receive array to M transmitted waveforms. Therefore, the effective aperture of 
the virtual array can be significantly extended that leads to improved angular resolution. In 
the methods of liTTl - lfTOl , the rotational invariance property is achieved by partitioning the 
receiving array to two (in ESPRIT case) or multiple (in PARAFAC case) overlapped subarrays. 
Then, the rotational invariance is also presumed for the virtual enlarged array of MN virtual 
antennas. However, the methods of [fT7l - [fT9l employ full waveform diversity, i.e., the number of 
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transmitted waveforms equals the number of transmit antennas, at the price of reduced transmit 
energy per waveform. In other words, for fixed total transmit energy, denoted hereafter as E, each 
waveform has energy E/M transmitted omni-directionally. It results in a reduced signal-to-noise 
ratio (SNR) per virtual antenna. On the other hand, it is well-known that the estimation accuracy 
of subspace-based techniques suffers from high SNR threshold, that is, the root-mean-squared 
error (RMSE) of DOA estimates approaches the corresponding Cramer-Rao bound (CRB) only 
for relatively high SNRs [|20l, Jill. Therefore, the higher the SNR, the better the DOA estimation 
performance can be achieved. 

It has been shown in [[22|. [|23l that a tradeoff between the SNR gain and aperture of the 
MIMO radar virtual array can be achieved by transmitting less than M orthogonal waveforms. 
Exploiting this tradeoff, we develop in this papeo a group of DOA estimation methods, which 
allow to increase the SNR per virtual antenna by (i) transmitting less waveforms of higher energy 
and/or (ii) focusing transmitted energy within spatial sectors where the targets are likely to be 
located. At the same time, reducing the number of transmitted waveform reduces the aperture 
of the virtual array, while a larger aperture may be useful for increasing the angular resolution 
at high SNR region. 

The contributions of this work are based on the observation that by using less waveforms 
the energy available for each transmitted waveform can be increased, that is, the SNR per each 
virtual antenna can be improved, while the aperture of the virtual array decreases to KN, where 
i^ < M is the number of orthogonal waveforms. Moreover, the SNR per virtual antenna can 
be further increased by focusing the transmitted energy in a certain sector where the targets 
are located. Our particular contributions are as follows. The transmit beamformers are designed 
so that the transmitted energy can be focussed in a certain special sector where the targets 
are likely to be located that helps to improve the SNR gain at each receive antenna, and 
therefore, improve the angular resolution of DOA estimation techniques, such as, for example, 
MUSIC -based techniques. Moreover, we consider the possibility of obtaining the rotational 
invariance property while transmitting 1 < K < M orthogonal waveforms using different 
transmit beamforming weight vectors and focusing the transmitted energy on a certain spatial 
sector where the targets are located. It enables us to design search-free ESPRIT/R^RAFAC -based 
DOA estimation techniques. In addition, we derive CRB for the considered DOA estimation 

'An early exposition of a part of this work has been presented in 1241 . 
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schemes that aims at further demonstrating how the DOA estimation performance depends on 
the number of transmitted waveforms, transmit energy focusing, and effective array aperture of 
the MIMO radar virtual array. 

The paper is organized as follows. In Section II, MIMO radar signal model is briefly in- 
troduced. In Section III, we present the transmit beamspace-based MIMO radar signal model. 
Two approaches for designing the transmit beamspace weight matrix are given in Section IV. 
In Section V, we present MUSIC and ESPRIT DOA estimation for transmit beamspace-based 
MIMO radar. Performance analysis and CRB are given in Section VI. Simulation results which 
show the advantages of the proposed transmit beamspace-based MIMO radar DOA estimation 
techniques are reported in Section VII followed by conclusions drawn in Section VIII. 

II. MIMO Radar Signal Model 

Consider a MIMO radar system equipped with a transmit array of M colocated antennas and 
a receive array of A^ colocated antennas. Both the transmit and receive antennas are assumed to 
be omni-directional. The M transmit antennas are used to transmit M orthogonal waveforms. 
The complex envelope of the signal transmitted by the mth transmit antenna is modeled as 

SnXt) = y^Mt), m = l,...,M (1) 

where t is the fast time index, i.e., the time index within one radar pulse, E is the total 
transmitted energy within one radar pulse, and 0m(i) is the mth baseband waveform. Assume 
that the waveforms emitted by different transmit antennas are orthogonal. Also, the waveforms 
are normalized to have unit-energy, i.e., /^ |0m(t)p(it = 1, m = 1, . . . , M, where T is the pulse 
width. 

Assuming that L targets are present, the A^ x 1 received complex vector of the receive array 
observations can be written as 

L 

x(t,r) = J]n(t,r)b(^z) + z(t,r) (2) 

1=1 

where r is the slow time index, i.e., the pulse number, h{9) is the steering vector of the receive 
array, z(t, r) is A^ x 1 zero-mean white Gaussian noise term, and 




is the radar return due to the Ith target. In ([3]), ai(r), 9i, and a(6';) are the reflection coefficient 
with variance a^, spatial angle, and steering vector of the transmit array associated with the /th 
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target, respectively, (j)(t) = [0i(t), . . . ,0A/(t)]^ is the waveform vector, and (■)^ stands for the 
transpose. Note that the reflection coefficient a/(r) for each target is assumed to be constant 
during the whole pulse, but varies independently from pulse to pulse, i.e., it obeys the Swerling 
II target model fil9i1 . 

Exploiting the orthogonality property of the transmitted waveforms, the A^ x 1 component of 
the received data Q due to the mth waveform can be extracted using matched- filtering which 
is given as follows 

x^(r) ^ J x(t, T)(l)*^{t)dt, m = 1, . . . , M. (4) 

where (■)* is the conjugation operator. Stacking the individual vector components dH) in one 
column vector, we obtain the MN x 1 virtual data vector ^ 

yMiMo(r) = [xf (r) ■ ■ ■ ^M{T)f 
L 



yjX]«z(T)a(^0®b(^0 + z(r) 




1=1 

— ^ai{T)uuiMo{Oi) + z{T) (5) 

1=1 

where ® denotes the Kronker product, 

UMiMoW = a(^)®b(0) (6) 

is the MN x 1 steering vector of the virtual array, and z(r) is the MN x 1 noise term whose co- 
variance is given by cx^^Imn- The MNxMN covariance matrix Rmimo — E{yMiMo(T")yMiMO (''")} 
is hard to obtain in practice. Therefore, the following sample covariance matrix 

1 '^ 

R-MIMO = TT / ,yMIMo(T)yMIMo(''") C^) 

^ r=l 

is used, where Q is the number of snapshots. 

III. Transmit Beamspace Based MIMO Radar Signal Model 

Instead of transmitting omin-directionally, we propose to focus the transmitted energy within 
a sector by forming K directional beams where an independent waveform is transmitted over 
each beam. Note that the spatial sector can be estimated in a preprocessing step using any 
low-resolution DOA estimation technique of low complexity. 

Let C = [ci, . . . , ck]'^ be the transmit beamspace matrix of dimension M x K (K < M), 
where c^ is the M x 1 unit-norm weight vector used to form the /cth beam. The beamspace matrix 
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can be properly designed to maintain constant beampattern within the sector of interest and to 
minimize the energy transmitted in the out-of-sector areas. Let 4>xi^) — [0i(^)) ■ ■ ■ ? ^K{t)V be 
the K X 1 waveform vector. The A;th column of C is used to form a transmit beam for radiating 
the kth waveform (pkii)- The signal radiated towards a hypothetical target located at a direction 
9 via the kth beam can be modeled as 

Skit,e) = ^{cfa{e))Mt) (8) 



where ^JE jK is a normalization factor used to satisfy the constraint that the total transmit 
energy is fixed to E. The signal radiated via all beams towards the direction Q can be modeled 
as 

fc=i 

= y|c^a(^)0^(t) = ^^^\Q)4>At) (9) 

where (•)^ stands for the Hermitian transpose and 

a(^)^(C^a(e))^. (10) 

Then, the transmit beamspace can be viewed as a transformation that results in changing the 
M X 1 transmit array manifold a(6') into the K x 1 manifold (flOl) . 

At the receive array, the A^ x 1 complex vector of array observations can be expressed as 

Xbeam(t,r) = J^^aK^) [^ {Ql)4>K{t)) h{d,)^'L{t,r) . (11) 

l=\ 

By matched-filtering Xbcam(^, t) to each of the waveforms </)fc (A; = 1, . . . , K), the received signal 
component associated with each of the transmitted waveforms can be obtained as 

yA;(r) = / Xbeam(t,^)0fc(^)c?i 

= A/f E«,(r)a[,](0,)b(0,) + z,.(r) (12) 

i=\ 

where [■)\k\ is the /cth entry of a vector and the A^ x 1 noise term is defined as 

Zk{T) = J z{t,T)cj)l{t)dt. (13) 
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Stacking the individual vector components ([72)) in one column vector, we obtain the following 
KN X 1 virtual data vector 



= V f E «'(^) (^(^') ® b(^0) + ~Mr) (14) 

1=1 

where zx(t) = [zi(r), . . . ,zx(t)]^ is the KN x 1 noise term whose covariance is given by 

2t 

The transmit beamspace signal model given by (fT4)) provides the basis for optimizing a general- 
shape transmit beampattern over the transmit beamspace weight matrix C. By carefully designing 
C, the transmitted energy can be focussed in a certain spatial sector, or divided between several 
disjoint sectors in space. As compared to traditional MIMO radar, the benefit of using transmit 
energy focusing is the possibility to increase in the signal power at each virtual array element. 
This increase in signal power is attributed to two factors: 

(i) transmit beamforming gain, i.e., the signal power associated with the /cth waveform reflected 
from a target at direction 6 is magnified by factor \c^ai{6)\'^\ 

(ii) the signal power associated with the A;th waveform is magnified by factor E/K due to 
dividing the fixed total transmit power E over K < M waveforms instead of M waveforms. 

As it is shown in subsequent sections, the aforementioned two factors result in increasing the 
SNR per virtual element which in turn results in lowering the CRB and improving the DOA 
estimation accuracy. 

IV. Transmit Beamspace Design 

In this section, two methods for designing the transmit beamspace weight matrix C are 
proposed. 

A. Motivations 

Given the beamspace angular sector-of-interest ©, several beamspace dimension reduction 
techniques applied to the data at the output of a passive receive array are reported in the literature 
(see ll25l and references therein). The essence of these beamspace dimension reduction techniques 
is to perform DOA estimation in the reduced dimension space rather than in the elementspace 
(full dimension of received data) which leads to great computational savings. Performing DOA 
estimation in a reduced dimension beamspace has also proved to improve probability of source 
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resolution as well as estimation sensitivity [|26ll . However, it is well known that the CRB on DOA 
estimation performed in reduced dimension space is higher than (or at best equal to) the CRB 
on DOA estimation performed in elementspace, i.e., applying beamspace dimension reduction 
techniques to passive arrays does not improve the best achievable estimation performance [27]. 
This is natural because beamspace dimension reduction in passive arrays just preserves the signals 
observed within a certain spatial sector and filters out signals observed outside that sector. 

The transmit beamspace processing proposed in Section Ull] focuses all transmit energy within 
the desired spatial sector instead of spreading it omni-directionally in the whole spatial domain. 
In other words, the amount of energy that otherwise could be wasted in the out-of- sector areas is 
added to the amount of energy to be transmitted within the desired sector. As a result, the signal 
strength within the desired sector is increased potentially leading to improved best achievable 
estimation performance, i.e., as will be seen in Section IVI-BI lowering the CRB. The transmit 
beamspace weight matrix C can be designed such that the following two main requirement 
are satisfied: (i) Spatial distribution of energy transmitted within the desired sector is uniform; 
(ii) the amount of energy that is inevitably transmitted in the out-of-sector area is minimized. 
Here we present two methods for satisfying these two requirements. 

B. Spheroidal Sequences Based Transmit Beamspace Design 

Discrete prolate spheroidal sequences (DPSS) have been proposed for beamspace dimension 
reduction in array processing ll28l . The essence of the DPSS-based approach to beamspace 
dimension reduction f2\\, J28l| is to maximize the ratio of the beamspace energy that comes 
from within the desired sector to the total beamspace energy, i.e., the energy within the 
whole spatial domain [—n/2, n/2]. Following this principle, we propose to design the transmit 
beamspace weight matrix so that the ratio of the energy radiated within the desired spatial sector 
to the total radiated energy is maximized. That is, the transmit beamspace matrix C is designed 
based on maximizing 



Jt J-f 


cfa 


{0)MtW 


dOdt 
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where A = J^ a{9)a^{9)d9 is a non-negative matrix. Note that the first equality in (fT5l) follows 
from the fact that J^ (j)(t)(j)* {t)dt = 1 . Moreover, it can be shown that if a(^) obeys Vandermonde 
structure, then the following holds ETTl 

/E. 
' |cfa(^)|2rf^ = 27rcfcfc. (16) 

-f 
Substituting (fT6l) in (fT5l) . we obtain 

r. = ^^, k = i,...,K. (17) 

The maximization of the ration in (flTl) is equivalent to maximizing its numerator while fixing its 
denominator by imposing the constraint ||cfc|| = 1. To avoid the trivial solution Ci = . . . = c^, 
an additional constraint might be necessary, for example, the orthogonality constraint c^^^Cfc/ = 
0, k ^ k' can be imposed. Then, the maximization of {Tk}^=i subject to the constraint C^C = I 
corresponds to finding the eigenvectors of A that are associated with the K largest eigenvalues 
of A. That is, the transmit beamspace matrix is given as 

C = [ui,U2,...,Ui^] (18) 

where {ui}f^^ are K principal eigenvectors of A. 

It is worth noting that any steering vector a(9),9 E & belongs to the space spanned by the 
columns of C. This means that the magnitude of the projector of a(9) onto the column space 
of C is approximately constant, i.e., the transmit power distribution 

H{9) = a"{9)CC"a (19) 

is approximately constant V^ G ©. Therefore, the total transmit power distribution, i.e., the 
distribution of the power summed over all waveforms, within the spatial sector is approxi- 
mately uniform. However, the transmit power distribution over individual waveforms may not be 
uniform within 0. Such a uniform transmit power distribution over individual waveforms may 
be desired especially when the difference between the first and the last essential eigenvalues 
of A is significant. In this respect, note that the orthogonality constraint C^C = I is imposed 
above only to avoid the trivial solution Ci = . . . = ck when maximizing (fTTI) . This is different 
from the case of the traditional beamspace dimension reduction techniques where orthogonality 
is required to preserve the white noise property. Therefore, in our case where the latter is not 
an issue, C can be easily rotated, if necessary, so that it achieves desirable features such as 
uniform transmit power distribution over individual waveforms. The orthogonality may be lost 
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after such rotation of C without any consequences. For example, such rotation of C is used in 
our simulations (Section VII) to ensure that the uniform transmit power distribution per each 
individual waveform is achieved. 

C. Convex Optimization Based Transmit Beamspace Design 

A more general transmit beamspace design technique applicable to a transmit array of arbitrary 
geometry can also be formulated. Moreover, some DOA estimation techniques can be applied 
only if the received data enjoy certain properties. For example, the ESPRIT DOA estimation 
technique requires that the received data snapshots consist of two sets where one data set is 
equivalent to the other up to some phase rotation. This property is commonly referred to as 
rotational invariance. More generally, PARAFAC-based DOA estimation techniques are based 
on the assumption that multiple received data sets are related via rotational invariance properties. 
Fortunately, the transmit beamspace matrix C can be designed so that the virtual data vectors 
given by (fT2)) enjoy such rotational invariance property. 

To obtain the rotational invariance property, we propose a convex optimization based method 
for designing C. For the virtual snapshots {yA;(T)}^i to enjoy rotational invariance, the following 
relationship should be satisfied 

em(e)|^cfa(^))|b(^) = e^^fc'(^)|(c^a(^))|b(0), V^ G 0, k,k'=l,...,K (20) 

where j = y/—T, and fXk{0) and fik'id) are arbitrary phase shifts associated with the fcth and k'th 
virtual snapshots, respectively. The relationship (|20|) can be satisfied if the transmit beamspace 
matrix C is designed/optimized to meet the following requirement 

C^a(e) = d(^), V0 e &. (21) 

where d{6) = [e^^'^^^\ . . . ,e^^^^^^]'^ is of dimension K x 1. A meaningful formulation of a 
corresponding optimization problem is to minimize the norm of the difference between the left- 
and right-hand sides of (|2TI) while keeping the worst transmit power distribution in the out-of- 
sector area below a certain level. This can be mathematically expressed as follows 

minmax||C^a(0i)-d(0)||, 9, e e, i = 1, . . . , I (22) 

C i 

subject to ||C^a(%)|| < 7, 9j e @, j = 1, . . . , J (23) 

where combines a continuum of all out-of-sector directions, i.e., directions lying outside the 
sector-of- interest 0, and 7 > is the parameter of the user choice that characterizes the worst 
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acceptable level of transmit power radiation in the out-of-sector region. The parameter 7 has an 
analogy to the stop-band attenuation parameters in the classic bandpass filter design and can be 
chosen in a similar fashion ||29l . 

V. Transmit Beamspace Based DOA Estimation 

In this section, we design DOA estimation methods based on transmit beamspace processing 
in MIMO radar. We focus on subspace based DOA estimation techniques such as MUSIC and 
ESPRIT. 

A. Transmit Beamspace Based MUSIC 

The virtual data model (O can be rewritten as 

ybcam(r) = Va(r) + zx(r) (24) 

where 

a(r)^ [«i(r),...,ai(r)]^ (25) 

V^[v(ei),...,v(ei)] (26) 

vW^y|(C^a(e))®b(e). (27) 
The KN X KN transmit beamspace-based covariance matrix is given by 

Rbeam = E {ybcam(r)y^eam(^) } 

= VSV^ + allKN (28) 

where S = E{Q:(r)a^(r)} is the covariance matrix of the reflection coefficients vector. The 
sample estimate of (|28] ) takes the following form 

1 "^ 

Rbeam = T^Yl ybeam(r)y^eam(^)- (29) 

^ r=l 



The eigendecomposition of (|29l ) can be written as 

R-beam = EgAgEg + EnAnE^ (30) 

where the L x L diagonal matrix As contains the largest (signal-subspace) eigenvalues and the 
columns of the KN x L matrix Eg are the corresponding eigenvectors. Similarly, the {KN — 
L) X [KN — L) diagonal matrix An contains the smallest (noise-subspace) eigenvalues while 
the KN X [KN — L) matrix En is built from the corresponding eigenvectors. 
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Applying the principle of the elementspace MUSIC estimator [|5]|, the transmit beamspace 
spectral-MUSIC estimator can be expressed as 

^^^) - V^(0)QV(0) ^^^^ 

where Q = EnE^ = I — EgEf^ is the projection matrix onto the noise subspace. Substituting 
(l27t into dSB, we obtain 

^ [(C^a(g))^b(g)]^[(C^a(g))^b(g)] 

[(C^a(^)) ® h{e)f Q [(C^a(^)) ® b 
_ [a^(g)CC^a(g)] ■ [b^(g)b(g)] 
~ [(C^a(0)) ® b(^)]^ Q [(C^a(^)) ® b 
^ iVa^(g)CC^a(g) 

[(C^a(0)) b(^)]^ Q [(C^a(^)) ® b(0)] ' 

5. Transmit Beamspace Based ESPRIT 

The signal component of the data vectors {yk{T)}^^i can be expressed as 

yfc(r) = T,a(r) (33) 

where 

T, ^[b(^i),...,b(^i)]*fc (34) 

*fc ^ diag {cf a(^i), . . . , cf a(0i)} . (35) 

It is worth noting that the matrices {Tfc}^^ are related to each other as 

Tfc = T,*7i*fc, k,j = l,...,K. (36) 

By carefully designing the beamspace weight matrix C, for example by using (|22)) - (l23l) . the 
relationship (|36l) enjoys the rotational invariance property. 

Consider the case when only two transmit beams are formed. Then, the transmit beamspace 
matrix is C = [ci,C2]. In this case, (|36l) simplifies to 

T2 = Ti* (37) 

where 



Furthermore, ( |38] ) can be rewritten as 

* = diag{A{9i)e'^^^'\...,A{eL)e'^^^^^} (39) 
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where A(9) and n(9) are the magnitude and angle of c^a(6')/cf^a(6'), respectively. It can be 
seen from (l37l) and (|39| ) that the vectors yi and y2 (see also (1331) ) enjoy the rotational invariance 
property. Therefore, ESPRIT-based DOA estimation techniques can be used to estimate ^ and 
{9i}JLi can be obtained from ^ by looking up a table that converts Vl{9) to 9. Moreover, if 
i^ > 2 is chosen, more than two data sets which enjoy the rotational invariance property can be 
obtained by properly designing the transmit beamspace weight matrix. In this case, the means 
of using PARAFAC instead of ESPRIT are provided as well. 

It is worth mentioning that for traditional MIMO radar ([5]), DOA estimation using ESPRIT 
has been proposed in ifTTl . Specifically, y(r) in ^ has been partitioned into yi(r) = [xf (r), 
. . . ,xj^_]^(r)]-^ and y2(T) = [x|'(r), . . . ,xfj(r)]-^ and it has been shown that yi and y2 
obey the rotational invariance property which enables the use of ESPRIT for DOA estimation. 
However, the rotational invariance property is valid in [I17II only when the transmit array is 
a uniform linear array (ULA). Therefore, the method of ifTTl is limited by the transmit array 
structure and may suffer from performance degradation in the presence of array perturbation 
errors. Moreover, it suffers from low SNR per virtual antenna as a result of dividing the total 
transmit energy over M different waveforms. 

VI. Performance Analysis 

In this section, we analyze the performance of the proposed transmit beamspace MIMO radar 
DOA estimation approach as compared to the MIMO radar DOA estimation technique. For 
the reason of comparison, we also consider two techniques that have been recently reported 
in the literature that employ the idea of dividing the transmit array into several smaller sub- 
apertures/subarrays. The first technique uses transmit subapertures (TS) for omni-directionally 
radiating independent waveforms [|30ll . [|3T1l . If the number of subapertures is chosen as K < M, 
then each subaperture radiates pulses of energy E/K. The second technique is based on par- 
titioning the transmit array into overlapped subarrays where the antennas that belong to each 
subarray are used to coherently transmit an independent waveform [|22ll . [|23l . We refer to this 
technique as transmit array partitioning (TAP). Note that the TAP technique has transmit coherent 
gain while the TS technique does not have such a coherent transmit gain. 

In the following subsection, we compare between the aforementioned techniques in terms 
of the effective aperture of the corresponding virtual array, the SNR gain per virtual element, 
and the computational complexity associated with eigendecomposition based DOA estimation 

July 5, 2010 DRAFT 



14 

techniques. In the next subsection, we express/discuss the CRB for all considered techniques. 

A. Colocated Uniform Linear Transmit Array 

Consider the case of a ULA at the transmitter with A/2 spacing between adjacent antennas, 
where A is the propagation wavelength. Taking the first antenna as a reference, the transmit 
steering vector can be expressed as 

a(0) 4 [1, e-J''^^'"^, . . . , (,-MM-i)sineY ^ ^^^^ 

Also, the receive antennas are assumed to be grouped in a ULA with half a wavelength interele- 
ment spacing. Then, the receive steering vector is given as 

b(^)4 [l,e-^"^^'^^...,e-^"(^^i)^'"^]^. (41) 

It is worth noting that the transmit and receive array apertures are (M — l)A/2 and (A^ — l)A/2, 
respectively. In light of (l40l) and (|4T]) . we discuss/analyze the following cases. 

1) Traditional MIMO radar: Substituting (l40l) and (|4T1) in ^, the MIMO radar virtual steering 
vector can be expressed as 

j-K{N—l)s\r\8 —JTTsinO —jnN sin9 



UMiMoW = [l,...,e^^'^(^-l)^-^e 



. . .,e 



^_jV(M-l)sine _ _ _ g-j27r(A/-l+Af-l)sinei^ _ ^42) 

From (|42|) . we observe that the effective virtual array aperture is (M + A^ — 2) A/2. Note that 
the virtual steering vector a(9) (g) a{9) is of dimension MN x 1, yet it only contains M + N — 1 
distinct elements. Moreover, the SNR gain per virtual element is proportional in this case to 
E/M. This low SNR gain can lead to poor DOA estimation performance especially at low SNR 
region. The computational complexity of applying eigen-decomposition based DOA estimation 
techniques is of 0{M^N^) in this case. 

2) Transmit subaperturing based MIMO radar: As compared to the traditional MEMO radar, 
TS-based MIMO radar employs K subapertures instead of M. This results in higher SNR per 
virtual element in the corresponding virtual array at the receiver. To capitalize on the effect of 
this factor on the DOA estimation performance, we consider the extreme case when only two 
transmit subaperture^ are used to radiate the total transmit energy E. In this case, each transmit 

One can think of a subaperture as a large omni-directional antenna which is capable of radiating energy E /2 instead of 
E/M. This case might be practically unattractive as it will require power amplifier of much higher amplifying gain as compared 
to the case of using M transmit antennas. However, for the sake of theoretical analysis/comparison we consider it in this paper. 
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subaperture radiates a waveform of energy E /2. Assume that the two transmit subapertures are 
separated in space by C, wavelength and that the first transmit subaperture is taken as a reference. 
Then, the MIMO radar data model © becomes of dimension 2N x 1 and can be expressed as 

yTs(r) = a/I E«/(r) ([1, e-^-2<^-^f ®b(^)) +ZTs(r) (43) 

1=1 

where ZTs(r) = [zi(r), Z2(r)]^. 

Comparing ([5]) to (|43l) . we observe that the signal strength for MIMO radar with M transmit 



antennas is proportional to -^ EjM while the signal strength for the TS-based MIMO radar is 
proportional to \/E/2. This means that (gS]) offers an SNR gain that is M/2 times the SNR 
gain offered by ©. We also observe from (|43|) that the effective virtual array aperture is given 
by (C + ^ ~" l)'^/2- Therefore, the effective aperture for this case can be controlled by selecting 
C,. For example, selecting ( = \/2 yields a virtual array steering vector of dimension 2N that 
contains only A^ + 1 distinct elements, i.e, the effective aperture would be N\/2. This case 
is particularly important when performing DOA estimation using search free techniques such 
as ESPRIT. Another important case is the choice ( = N\/2 which yields a virtual array that 
is equivalent to a 2A^-element ULA. In this case, the effective array aperture will be (2A^ — 
l)A/2. The computational complexity of applying eigen-decomposition based DOA estimation 
techniques is of 0(2^A^'^) in this case. 

3) Transmit array partitioning based MIMO radar: Following the guidelines of [23] and 
selecting K = 2, the M-antenna transmit array is assumed to be partitioned into two fully 
overlapped subarrays of M — 1 antenna each. The (M — 1) x 1 beamforming weight vectors 
wi = W2 = w are used to form transmit beams that cover the spatial sector 0. Two independent 
waveforms are radiated. Each waveform has E/2 energy per pulse. The transmit weight vector 
w can be designed such that the transmit gain is approximately the same within the desired 
sector 0, i.e., |w^a(^)| = Gtap, V6' G where a{9) contains the first M — 1 elements of 
the vector a(9) and Gtap is the TAP transmit coherent processing gain. Then, the TAP-based 
MIMO radar data model becomes of dimension 2A^ x 1 and can be formally expressed as 

yTAp(r) = Y^X^"K^) (w^a(6'z)) utap(6'z) + ZTAp(r) 
1=1 

= GtapJ i^^ai{T)uTAp{6i) + ztap(t) (44) 

1=1 
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where utap(^) — [l, e^-''"''™^]^ ® h{9) is the 2N x 1 steering vector of the corresponding 
virtual array. Note that the TAP-based MIMO radar has transmit coherent gain Gtap which 
results in improvement in SNR per virtual element. However, the corresponding virtual array 
contains A^ + 1 distinct elements, i.e., the effective virtual array aperture is limited to NX/ 2. The 
computational complexity of applying eigen-decomposition based DOA estimation techniques is 
of 0{2^N^) in this case. 

4) Transmit beamspace based MIMO radar: For the proposed transmit beamspace MIMO 
radar, we choose K = 2 and use (|22|)-(|231) for designing C = [ci C2] such that 

'-'"beam ' y^) ^ J 7 



d(^) 



V^gG 



(45) 



where Gbcam is the transmit beamspace gain, i.e., |cf^a(0)| ^ |c^a(0)| ?^ Gbcam, V6' G 0. 
This yields a virtual array with 2N distinct elements and {2N — l)A/2 effective array aperture. 
Moreover, the proposed transmit beamspace technique offers SNR gain of Gbcam ■ E/M, i.e., it 
combines all the benefits of all other aforementioned techniques. The computational complexity 
of applying eigen-decomposition based DOA estimation techniques is of 0{2^N^) in this case. 
A comparison between all methods considered is summarized in the following table. 



Table 1: Comparison between transmit beamspace-based MIMO radar and other existing tech- 
niques. 





Effective aperture 


SNR gain per virtual element 


Computational complexity 


Traditional MIMO (5) 


(M + N -2)^ 


E 
M 


0{M^N'') 


Transmit subaperturing {C, = ^) 


Ar| 


B 
2 


0{2^N^) 


Transmit subaperturing {C, = A''-|) 


(2iV-l)A 


B 
2 


0{2^N^) 


Transmit array partitioning {AA\ 


N^ 


r^2 E 
'-'TAP ■ 2" 


0{2^N^) 


Transmit beamspace MIMO \45\ 


(2iV-l)| 


f^2 E 


0{2^N^) 



B. Cramer-Rao Bound 

In this section, we discuss the CRB on DOA estimation accuracy in transmit beamspace-based 
MIMO radar. 

In the case of transmit beamspace-based MIMO radar, the virtual data model (fT4l) satisfies 
the following statistical model: 



ybeam(r) ~ Mc (fJ'ir) , R) 



(46) 
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where He denotes the complex muhivariate circularly Gaussian probability density function, 
/^(r) is the mean of ybcam(T), and R is its covariance matrix. 

1) Stochastic CRB: The stochastic CRB on estimating the DOAs using the data model (fT4l) 
is derived by assuming /j,(r) = and R = Rbcam, where Rbeam is given by (|28l) . Note that 
under these assumptions, the virtual array signal model (fT4l) (or its equivalent representation 
(|24)) ) has the same form as the signal model used in ll32l to derive the stochastic CRB for DOA 
estimation in conventional array processing. Therefore, the CRB expressions in its general form 
derived in [|32l can be used for computing the stochastic CRB for estimating the DOAs based 
on dill), that is. 



a? 



CRB(0) = ^ {Re(DP^D) G^} ' (47) 

where P^ — V(V^V)^^V^ is the projection matrix onto the space spanned by the columns 
of V, and G = (SV^R-^VS). In ^, D = [d(^i), . . . , (1(0^)] is the matrix whose Ith column 
is given by the derivative of the /th column of V with respect to 9i, i.e.. 



d9 \ K dO 

= wj ((C^a'(^)) ® b(^) + (C^a(0)) ® h'{d)) (48) 

where a'(0) = da{e)/de and b'(^) = dh{e)/de. 

2) Deterministic CRB: The deterministic CRB is derived by assuming At(r) = \'a.(T) and 
R = a^IxN- Under these statistical assumptions, the virtual data model (fT4)) is similar to the 
general model used in [[33l to find the deterministic CRB. According to [[33l . the deterministic 
CRB expression can be obtained from (l47l) by replacing G with S, where S is the sample 
estimate of S. 

It is worth noting that the expression (l47l) can be used not only for computing the CRB for the 
proposed transmit beamspace technique but also for the other techniques summarized in Table 1 . 
Indeed, the following cases show how these techniques can be viewed as special cases of the 
proposed model (fT4l) . 

1) Choosing C = Ij\/, the transmit beamspace signal model (fT4l) simplifies to the traditional 
MIMO radar signal model ([5]). Therefore, the CRB for the traditional MIMO radar can be 
obtained by substituting C = Im in (ST]) and (|48|) . 

2) The TS-based MEMO radar signal model with (^ = A/2 can be obtained from (fT4l) by 
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choosing C in the following format 



(49) 



1 
1 

3) The TS-based MIMO radar signal modal for ( = N\/2 can be obtained from (fT4)) by 
choosing C in the following format 



(50) 



1 
1 
4) Finally, the TAP-based MIMO radar signal model (1441) can be obtained from (fT41) by 
choosing C in the following format 

wi 

(51) 

W2 

VII. Simulation Results 

Throughout our simulations, we assume a uniform linear transmit array of M = 10 omni- 
directional antennas spaced half a wavelength apart. At the receiver, A^ = 10 omni-directional 
antennas are also assumed. The additive noise is Gaussian zero-mean unit-variance spatially and 
temporally white. Two targets are located at directions —1° and 1°, respectively. The sector of 
interest = [—5°, 5°] is taken. Several examples are used to compare the performances of the 
following methods: (i) The traditional MIMO radar ©; (ii) The TS-based MIMO radar dUl) 
with C = A/2; (iii) The TS-based MIMO radar ^ with C = A^A/2; (iv) the TAP-based MIMO 
radar (|5TI) : and (v) The proposed transmit beamspace MIMO radar ([14)) . For all methods tested, 
the total transmit energy is fixed to E = M. For the traditional MIMO radar ([5]), each transmit 
antenna is used for omni-directional radiation of one of the baseband waveforms 

le^2.ft^ m = l,...,M. (52) 

For all methods that radiate two waveforms only, the first and second waveforms of (|52l) are 
used. The total number of virtual snapshots used to compute the sample covariance matrix is 
Q = 300. In all examples, the RMSEs and the probability of source resolution for all methods 
tested are computed based on 500 independent runs. 



(pm(t) 
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A. Example 1 : Stochatstic and Deterministic CRBs 

For the proposed spheroidal sequences based transmit beamspace MIMO radar (fTSi) . the non- 
negative matrix A = j^ai{0)ai^ {6)d6 is built. The two eigenvectors associated with the largest 
two eigenvalues are taken as the principle eigenvectors, i.e., C = [ui U2] is used. Two waveforms 
are assumed to be radiated where each column of C is used to form a transmit beam for radiating 
a single waveform. Fig. [T] shows the transmit power distribution of the individual waveforms 
|cfca(^)p, fc = 1, 2 as well as the distribution of the total transmitted power (fT9l ). As we can see 
from this figure, the individual waveform power is not uniformly distributed within the sector 
© while the distribution of the total transmitted power is uniform. To achieve uniform transmit 
power distribution for each waveform a simple modification that involves vector rotation to the 
transmit beamspace matrix C can be performed. This means that C can be modified as 

C = CQ (53) 

where Q is a 2 x 2 unitary matrix defined as 

Q = 



(54) 



VV2 VV2 

v/T72 -v/I72_ 

Fig. |2] shows the transmit power distribution for the individual waveforms |cfca(6')p, /c = 1, 2 as 
well as for the total transmitted power (fT9l) . It can be seen from this figure that the distribution 
of individual waveforms is uniform within the desired sector. 

For the TAP-based method, the transmit array is partitioned into two overlapped subarrays 
of 9 antennas each. Each subarray is used to focus the radiation of one waveform within the 
sector 0. The transmit weight vectors are chosen as wi = W2 = [—0.5623 —0.5076 —0.4358 
-0.3501 -0.2542 -0.1524 -0.0490 0.0512 0.1441]^. This specific selection is obtained by 
averaging the two eigenvectors associated with the maximum two eignevalues of the matrix 
Jq ai(^)af (6')(i6', where ai(0) is the 9x1 steering vector associated with the first subarrayO 
The transmit power distribution of both waveforms is exactly the same as shown in Fig. |3l 

The stochastic CRBs for all methods considered are plotted versus SNR = cr^/o"^ in Fig. |4l It 
can be seen from this figure that the TS-based MIMO radar with C, = A/2 has the highest/worst 
CRB as compared to all other methods. Its poor CRB performance is attributed to the omni- 
directional transmission, i.e., wasting a considerable fraction of the transmitted energy within 

^Note that if the transmit array is not a ULA, then wi and W2 can be designed independently using classic FIR filter design 
techniques. 
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the out-of-sector region, and to the small effective aperture of the corresponding virtual array. 
We can also see from the figure that the TS-based MIMO radar with (^ = N\/2 exhibits much 
lower stochastic CRB as compared to the case with ( = A/2. The reason for this improvement 
is the larger effective aperture of the corresponding virtual array. The traditional MIMO radar 
with M transmit antennas has the same effective aperture as that of the TS-based MIMO radar 
with ( = N\/2 but lower SNR per virtual element. Therefore, the stochastic CRB for traditional 
MIMO radar is higher than the CRB for the TS-based MIMO radar with C = A^A/2. At the 
same time, it is better than the CRB of the TS-based MIMO radar with ( = A/2 due to larger 
effective aperture. The TAP-based MIMO radar has the same effective aperture as that of the TS- 
based MIMO radar with (^ = A/2 but higher SNR per virtual element due to transmit coherent 
processing gain. It yields lower CRB. In fact, the CRB for the TAP-based MIMO radar is 
comparable to that of the traditional MIMO radar. Finally, the transmit beamspace MIMO radar 
with spheroidal sequences based transmit weight matrix has the lowest CRB as compared to all 
other methods. This can be attributed to the fact that the proposed transmit beamspace-based 
MIMO radar combines the benefits of having high SNR due to energy focusing, high power of 
individual waveforms, and large effective aperture of the corresponding virtual array. 

The deterministic CRBs for all methods considered are plotted in Fig. [5l As can be seen from 
this figure, the same observations and conclusion that are drawn from the stochastic CRB curves 
also apply to the deterministic CRB. 

B. Example 2: MUSIC-based DOA Estimation 

In this example, the MUSIC algorithm is used to estimate the DOA for all aforementioned 
methods. Note that the targets are considered to be resolved if there are at least two peaks in 
the MUSIC spectrum and the following is satisfied [|21|| 



ei-Oi 



<^, / = 1,2 (55) 



where A6' = \62 — Oi\. Fig. |6] shows the probability of source resolution versus SNR for all 
methods tested. It can be seen from this figure that all methods exhibit a 100% correct source 
resolution at high SNR values. As the SNR decreases, the probability of source resolution starts 
dropping for each method at a certain point until it eventually becomes zero. The SNR level 
at which this transition happens is known as SNR threshold. It can be seen from Fig. [6] that 
the TS-based MIMO radar with C, = A/2 has the highest SNR threshold while the traditional 
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MIMO radar and the TAP-based MIMO radar have the second and third highest SNR thresholds, 
respectively. The SNR threshold of the TS-based MIMO radar with (^ = N\/2 is lower than the 
aforementioned three methods while the proposed transmit beamspace-based MIMO radar has 
the lowest SNR threshold, i.e., the best probability of source resolution performance. 

Fig. Ul shows the RMSEs for the MUSIC -based DOA estimators versus SNR for all methods 
tested. It can be seen from this figure that the TS-based MIMO radar with ( = A/2 has 
the highest/poorest RMSE performance. It can also be seen that the TAP-based MIMO radar 
outperforms the traditional MIMO radar at low SNR region while the opposite occurs at high 
SNR region. This means that the influence of having large effective aperture is prominent at 
high SNR region, while the benefit of having high SNR gain per virtual antenna (even if the 
effective aperture is small) is feasible at low SNR region. It can also be observed from Fig. |7] 
that the estimation performance of the TS-based MIMO radar with ( = N\/2 is better than 
that of both the traditional MIMO radar and the TAP-based MIMO radar. Finally, the proposed 
transmit beamspace-based MIMO radar outperforms all aforementioned methods. 

It is worth noting that the width of the desired sector is 10°. Therefore, the parts of the RMSE 
curves where the RMSEs exceed 10° in Fig. |7]are not important. Thus, the comparison between 
different methods within that region is meaningless. 

C. Example 3: ESPRIT-based DOA Estimation 

In this example, all parameters for all methods are the same as in the previous example except 
for the M X 2 transmit weight matrix associated with transmit beamspace-based MIMO radar 
which is designed using (|22|)-(|231). The out-of-sector region is taken as = [—90°, —15°] U 
[15°, 90°] and the parameter that controls the level of radiation within is taken as 7 = 0.38. 
Each column of the resulting matrix C is scaled such that it has unit norm. The transmit power 
distribution for this case is similar to the one shown in Fig. |2]and, therefore, is not shown here. 
ESPRIT-based DOA estimation is performed for all aforementioned methods. For the traditional 
MIMO radar-based method, the MN x 1 virtual array is partitioned into two overlapped subarrays 
of size (M — 1)A^ x 1 each, i.e., the first subarray contains the first (M — 1)A^ elements while 
the second subarray contains the last (M — 1)A^ elements. For all other methods, the 2N x 1 
virtual array is partitioned into two non-overlapped subarrays, i.e., the first subarray contains the 
first N elements while the second subarray contains the last A^ elements. 

The probability of source resolution and the DOA estimation RMSEs versus SNR are shown 
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for all methods tested in Figs.[8]and|9l respectively. It can be seen from Fig. [8] that the ESPRIT- 
based DOA estimator for the TS-based MIMO radar with ( = A/2 has the highest/poorest SNR 
threshold while the ESPRIT-based DOA estimator for the traditional MIMO radar has the second 
highest SNR threshold. The ESPRIT-based DOA estimator for the TAP-based MIMO radar has 
SNR threshold that is lower than the previous two estimators. Moreover, the SNR threshold of 
the ESPRIT-based DOA estimator for the TS-based MIMO radar with C = A^A/2 is lower than 
the SNR threshold of the ESPRIT-based DOA estimator for the TAP-based MIMO radar. Finally, 
the ESPRIT-based DOA estimator for the proposed transmit beamspace-based MIMO radar has 
the lowest SNR threshold, i.e., the best probability of source resolution performance. 

It can be seen from Fig. |9] that the ESPRIT-based DOA estimator for the TS-based MIMO 
radar with (^ = A/2 has the highest/poorest RMSE performance. Also this figure shows that the 
ESPRIT-based DOA estimator for the TAP-based MIMO radar outperforms the ESPRIT-based 
DOA estimator for the traditional MIMO radar at low SNR region while the opposite occurs at 
high SNR region. This confirms again the observation from the previous example that having 
large effective aperture is more important at high SNR region while having high SNR gain per 
virtual antenna is more important at low SNR region. It can also be observed from Fig. |9] that the 
ESPRIT-based DOA estimator for the TS-based MIMO radar with ( = NX/2 outperforms the 
ESPRIT-based DOA estimator for both the traditional MIMO radar and the TAP-based MIMO 
radar. Finally, the ESPRIT-based DOA estimator for the proposed transmit beamspace-based 
MIMO radar outperforms all aforementioned estimators. 

VIII. Conclusion 

A transmit beamspace energy focusing technique for MIMO radar with application to direction 
finding for multiple targets is proposed. Two methods for focusing the energy of multiple (two 
or more) transmitted orthogonal waveforms within a certain spatial sector are developed. The 
essence of the first method is to employ spheroidal sequences for designing transmit beamspace 
weight matrix so that the SNR gain at each receive antenna is maximized. The subspace 
decomposition-based techniques such as MUSIC can then be used for direction finding for 
multiple targets. The second method uses convex optimization to control the amount of dissi- 
pated energy in the out-of-sector at the transmitter and to achieve/maintain rotational invariance 
property at the receiver. This enables the application of search-free DOA estimation techniques 
such as ESPRIT. Performance analysis of the proposed transmit beamspace-based MIMO radar 
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and comparison to existing MIMO radar techniques with colocated antennas are given. Stochastic 
and deterministic CRB expressions as functions of the transmit beamspace weight matrix are 
found. It is shown that the proposed technique has the lowest CRB as compared to all other 
techniques. The computational complexity of the proposed method can be controlled by selecting 
the transmit beamspace dimension, i.e., by selecting the number of transmit beams. Simulation 
examples show the superiority of the proposed technique over the existing techniques. 

References 

[1] J. Li and P. Stoica, MIMO Radar Signal Processing. New Jersy: Wiley, 2009. 

[2] A. Haimovich, R. Blum, and L. Cimini, "MIMO radar with widely separated antennas," IEEE Signal Processing Magaz., 

vol. 25, pp. 116-129, Jan. 2008. 
[3] J. Li and R Stoica, "MIMO radar with colocated antennas," IEEE Signal Processing Magaz., vol. 24, pp. 106-114, Sept. 

2007. 
[4] V. R Pisarenko, "The retrieval of harmonics from a covariance function," Geophys. J. Royal Astronomical Soc, vol. 33, 

pp. 347-366, 1973. 
[5] R. O. Schmidt, "Multiple emitter location and signal parameter estimation," in Proc. RADC Spectral Estimation Workshop, 

Rome, NY, 1979, pp. 234-258. 
[6] R. O. Schmidt, "Multiple emitter location and signal parameter estimation," IEEE Trans. Antennas and Propagation, 

vol. 34, pp. 276-280, Mar. 1986. 
[7] A. J. Barabell, "Improving the resolution performance of eigenstructure-based direction-finding algorithms," in Proc. IEEE 

Int. Conf. on Acoustics, Speech and Signal Processing, Boston, MA, Apr. 1983, pp. 336-339. 
[8] Y. Bresler and A. Macovski, "Exact maximum likelihood parameter estimation of superimposed exponential signals in 

noise," IEEE Trans. Acoustics, Speech, and Signal Processing, vol. 34, no. 10, pp. 1081-1089, Oct. 1986. 
[9] R. Roy and T. Kailath, "ESPRIT - estimation of signal parameters via rotational invariance techniques," IEEE Trans. 
Acoustics, Speech, and Signal Processing, vol. 37, no. 7, pp. 984-995, July 1989. 
[10] R Li and R. J. Vaccaro, "Analysis of min-norm and MUSIC with arbitrary array geometry," IEEE Trans. Aerospace and 

Electronic Systems, vol. 26, no. 6, pp. 976-985, Nov. 1990. 
[II] M. Viberg and B. Ottersten, "Sensor array processing based on subspace fitting," IEEE Trans. Signal Processing, vol. 39, 

no. 5, pp. II 10-1121, May 1991. 
[12] M. Haardt, Efficient One-, Two-, and Multidimentional High-Resolution Array Signal Processing, Ph.D. Thesis, Shaker 

Verlag, Aachen, 1997. 
[13] N. D. Sidiropoulos, R. Bro, and G. B. Giannakis, "Parallel factor analysis in sensor array processing," IEEE Trans. Signal 

Processing, vol. 48, no. 8, pp. 2377-2388, Aug. 2000. 
[14] S. A. Vorobyov, Y. Rong, N. D. Sidiropoulos, and A. B. Gershman, "Robust iterative fitting of multilinear models," IEEE 

Trans. Signal Processing, vol. 53, no. 8, pp. 2678-2689, Aug. 2005. 
[15] L. Xu, J. Li, and P. Stoica, "Target detection and parameter estimation for MIMO radar systems," IEEE Trans. Aerosp. 

Electron. Syst., vol. 44, no. 3, pp. 927-939, July 2008. 
[16] I. Bekkerman and J. Tabrikian, "Target detection and localization using MIMO radars and sonars," IEEE Trans. Signal 
Processing, vol. 54, no. 10, pp. 3873-3883, Oct. 2006. 



July 5, 2010 DRAFT 



24 

[17] C. Duofang, C. Baixiao, and Q. Guodong, "Angle estimation using ESPRIT in MIMO radar," Electronics Letters, vol. 44, 

no. 12, pp. 11Q-11\, June 2008. 
[18] C. Jinli, G. Hong, and S. Weimin, "Angle estimation using ESPRIT without pairing in MIMO radar," Electronics Letters, 

vol. 44, no. 24, pp. 1422-1423, Nov. 2008. 
[19] D. Nion and N. D. Sidiropoulos, "A PARAFAC-based technique for detection and localization of multiple targets in a 

MIMO radar system," in Proc. IEEE Int. Conf. Acoustics, Speech and Signal Processing, Taipei, Taiwan, Apr. 2009, 

pp. 2077-2080. 
[20] P. Stoica and A. Nehorai, "MUSIC, maximum likelihood, and Cramer-Rao bound: further results and comparisons," IEEE 

Trans. Acoustics, Speech, and Signal Processing, vol. 38, no. 12, pp. 2140-2150, Dec. 1990. 
[21] H. L. Van Trees, Optimum Array Processing, Wiley, NY, 2002. 
[22] A. Hassanien and S. A. Vorobyov, "Transmit/receive beamforming for MIMO radar with colocated antennas," in Proc. 

IEEE Inter Conf. Acoustics, Speech, and Signal Processing, Taipei, Taiwan, Apr. 2009, pp. 2089-2092. 
[23] A. Hassanien and S. A. Vorobyov, "Phased-MIMO radar: A tradeoff between phased-array and MIMO radars," IEEE 

Trans. Signal Processing, vol. 58, no. 6, pp. 3137-3151, June 2010. 
[24] A. Hassanien and S. A. Vorobyov, "Direction finding for MIMO radar with colocated antennas using transmit beamspace 

preprocessing," in Proc. IEEE Inter. Workshop Computational Advances in Multi-Sensor Adaptive Processing, Aruba, Dutch 

Antilles, Dec. 2009, pp. 181-184. 
[25] A. Hassanien, S. Abd Elkader, A. B. Gershman, and K. M. Wong, "Convex optimization based beamspace preprocessing 

with an improved robustness against out-of-sector sources," IEEE Trans. Signal Processing, vol. 54, pp. 1587-1595, May 

2006. 
[26] H. B. Lee and M. S. Wengrovitz, "Resolution threshold of beamspace MUSIC for two closely spaced emitters," IEEE 

Trans. Acoust., Speech, Signal Processing, vol. 38, pp. 1545-1559, Sept. 1990. 
[27] S. Anderson, "On optimal dimension reduction for sensor array signal processing," Signal Processing, vol. 30, pp. 245-256, 

Jan. 1993. 
[28] P. Forster and G. Vezzosi, "Application of spheroidal sequences to array processing," in Proc. IEEE Inter Conf. Acoustics, 

Speech, and Signal Processing, Dallas, TX, May 1987, pp. 2268-2271. 
[29] B. A. Shenoi, Introduction to Digital Signal Processing and Filter Design, Wiley & Sons, Hoboken, New Jersey, 2006. 
[30] D. J. Rabideau and P. Parker, "Ubiquitous MIMO multifunction digital array radar," in Proc. 37th Asilomar Conf. on 

Signals, Systems, and Computer, Nov. 2003, vol. 1, pp. 1057-1064. 
[31] J. Bergin, S. McNeil, L. Fomundam, and P. Zulch, "MIMO phased-array for SMTI radar," in Proc. IEEE Aerospace Conf, 

Big Sky, MT, Mar. 2008, pp. 1-7. 
[32] P. Stoica, E. G. Larsson, and A. B. Gershman, "The stochastic CRB for array processing: A textbook derivation," IEEE 

Signal Processing Letters, vol. 8, pp. 148-150, May 2001. 
[33] P. Stoica and A. Nehorai, "Performance study of conditional and unconditional direction-of-arrival estimation," IEEE Trans. 

Acoust., Speech, Signal Processing, vol. 38, pp. 1783-1795 , Oct. 1990. 



July 5, 2010 DRAFT 



25 



20 



S 10 



CD 

cc 

H 

(n 

Q 

S -10 

5 

o 

Q_ 



-20- 



-30 



/ . I 



I I 
I I 
II 



- Transmit power distribution of(|)^(t) 
Transmit power distribution ofcl)^^ 

- Overall transmit power distribution 




I. I. . 1 .1 1 .h 



1,-1, 






, ,., '.I 1.1. 1. .1 .' I. . \ 



Ml 'M| 1 I 'i \ 

' I, ''n " ■', ^ 

I , '|. II : r :\ 

'i ■ '- i " ,1 . ^ Q 



-80 -60 -40 -20 20 40 60 80 

ANGLE (DEGREES) 



Fig. 1. Transmit beamspace beampattem using spheroidal sequences without rotation. Total transmit power is uniformly 
distributed within the desired spatial sector, however, different transmitted waveforms have different power distribution. 
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Fig. 2. Transmit beamspace beampattem using spheroidal sequences with rotation. Total transmit power as well as individual 
waveform powers are uniformly distributed within the desired spatial sector. 
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Fig. 3. TAP-based MIMO radar transmit beampattem using two fully overlapped subarrays. 
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Fig. 4. Stochastic CRB versus SNR. 
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Fig. 5. Deterministic CRB versus SNR. 
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Fig. 6. Probability of target resolution versus SNR for MUSIC-based DOA estimators. 
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Fig. 7. RMSE versus SNR for MUSIC-based DOA estimators. 
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Fig. 8. Probability of target resolution versus SNR for ESPRIT-based DOA estimators. 
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Fig. 9. RMSE versus SNR for ESPRIT-based DOA estimators. 
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